Impurity assisted nanoscale localization of plasmonic excitations in graphene 
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The plasmon modes of pristine and impurity doped graphene are calculated, using a real-space 
theory which determines the non-local dielectric response within the random phase approximation. 
A full diagonalization of the polarization operator is performed, allowing the extraction of all its 
poles. It is demonstrated how impurities induce the formation of localized modes which are absent in 
pristine graphene. The dependence of the spatial modulations over few lattice sites and frequencies 
of the localized plasmons on the electronic filling and impurity strength is discussed. Furthermore, it 
is shown that the chemical potential and impurity strength can be tuned to control target features of 
the localized modes. These predictions can be tested by scanning tunneling microscopy experiments. 
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Introduction: Graphene is an allotropic form of car- 
bon, in which atoms are arranged in a two-dimensional 
honeycomb lattice. Due to its thin crystal structure it is 
a promising material for a wide range of applications. [1| 
The prospects for its use in nanotechnology are fur- 
ther reinforced by the success of multiple fabrication 
alternatives. 0, Q Ab-initio calculations have shown that 
graphene is a gapless semiconductor, for which the va- 
lence and conduction bands meet at the Fermi energy. At 
this point the energy dispersion of the quasi-particle ex- 
citations has been found to be linear. Q Experimentally, 
the system can be tuned into a metallic regime by adjust- 
ing the chemical potential using a gate voltage. [5j Due to 
these unique properties there is considerable current in- 
terest in the electronic response of this material. 

While much theoretical research has focused on the 
ground state electronic properties of graphene, so far rel- 
atively few attempts have been made to study the col- 
lective excitations of this system. Nevertheless, there are 
some recent studies of the plasmonic modes in graphene 
using the linearized band structure [13, 
tight-binding band structure calculation, 
of a fully gapped band[l6| and doping 17 have also been 
explored. Some of the above mentioned papers have 
worked out the dispersion relation of plasmons. 

Recent STM experiments have shown that graphene 
is intrinsically disordered. fis| In order to determine its 
technological usefulness, it is therefore important to un- 
derstand the effects of disorder on its electronic proper- 
ties. Obviously, in the presence of impurities the system 
looses its translational symmetry, and it is not known 
what types of localized modes form around them, e.g. 
dipole, quadrupole, radial etc. Much recent research has 
focused on the effects of impurities on the ground state 
properties of graphene 1^-22, 2f| 2f|, but it is equally 
important to understand how impurities affect collective 
modes. In this letter, we examine the consequences of 
the presence of impurities on the plasmonic modes in 
graphene. Our main focus is to obtain and control lo- 
calized plasmons. The analysis of these features requires 
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the determination of their spatial profile, which - to the 
our best knowledge - has not yet been discussed in the 
literature on graphene. 

Model: The calculation of the dielectric response 
in metallic materials is conventionally performed using 
continuum- field Mie theory, [i?} However, such a semi- 
empirical continuum description breaks down beyond 
a certain de gree of roughness, introduced by atomic 
length scales. [28| Recently, a self-consistent quantum- 
mechanical approach was developed, which accounts for 
the non-locality of the dielectric response function. [29| 
Using this technique, the identification of plasmons was 
accomplished by scanning the frequencies of the modes 
with strongest induced fields. However, it does not pro- 
vide full information of all plasmon excitations supported 
by the system. In this contribution, we generalize this ap- 
proach. The polarization operator is diagonalized, pro- 
viding all its poles. Thus complete information of the 
plasmon excitations is obtained, including the local spec- 
tral densities of states. 

We model the electronic structure of graphene using a 
one-band tight-binding Hamiltonian, 
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where t = 2.7 'eV is the hopping parameter and n is the 
chemical potential. Uq is the magnitude of the impu- 
rity potential, x corresponds to its location, and a de- 
notes its spatial spread. U a is then the impurity po- 
tential felt at the other sites x a ^ xo, parametrized by 
U a = Uq exp ^ X °~ 2 X °^ j . In this article, the range of the 
impurity potential a will be only few lattice sites. H is 
diagonalized numerically, providing the eigenstates 
and eigenvalues E®. 

The direct Coulomb interaction is considered as a per- 
turbation, 
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where y>j(x) is the p x orbital at site j. In the lattice basis 
the induced charge is given by 
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Here the overlap of p z orbitals at different sites is ne- 
glected <£*(x)<£(,(x) = 5 a b(pl(x)(pb(x). The linear re- 
sponse charge equation in this basis is 



<WM = U ^nn (W) 
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Within the random phase approximation (RPA), the 
polarization operator is then obtained via H(lo) = 
n° (lo) (1 - Vn° (lo)) where n° (w) is the polarization 
operator of the non-interacting system. In the basis of 
eigenstates 1^"), n (lo) can be written as 
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This fourth rank tensor can be regarded as a matrix act- 
ing on the vectors tfi a b, and thus can be represented as 
the product 



n° (lo) = An (AE — wiy 1 , 
where An and AE are diagonal matrices 
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The polarization of the interacting system can then be 
expressed as 



n (lo) = An {AE -ul- VAn)" 1 
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Hence the plasmons correspond to charge densities bp 
An<f) such that 
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When the matrix AE — AnV is diagonalized, the po- 
larization has poles at lo = Ab for each eigenvalue At, of 
AE - AnV. 

The retarded Green's function H R (lo) = H(lo + i0 + ) 
gives the spectral density function A (lo) = 
— ^ImIl R (uj), which is expressed as 
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Due to the delta function A a p n $ (lo) is non-zero only 
for plasmonic frequencies. The only spatial profiles it 



<t>mn satisfying {AE — loI — V An) = 0. In order to 
plot the plasmon density of states the representation 
8 (lo — uj a ) — > J yi +1 2 ^ 1S use d, where 7 is chosen to 

be 0.05eV. Then the plasmon density of states trA(uo) 
can be written as 



trA (lo) = ^A a (uj) = ^A a - 
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A a will be referred to as the strength of the mode 
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where Anc/) a — Sp is the charge profile of the mode a. 
Using Eq. [121 the plasmon is considered to be localized 
when the major contribution to the sum comes from sites 
around the impurity. The following discussion of our re- 
sults is focussed on such localized modes. The plasmon 
density of states trA (lo) can be accessed experimentally 
through scanning tunneling experiments that reveal ei- 
ther direct or inelastic tunneling signatures associated 
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with plasmons or related lifetime effects. 

For the numerical simulation a finite-size realization of 
the graphene lattice with 96 sites is considered. An im- 
purity affects approximately one hexagon of the honey- 
comb lattice, as shown in Fig. QTa). Periodic boundary 
conditions are used to eliminate the boundary modes. 
LAPACK routines are used for the numerical diagonal- 
ization of Hq and n (lo). (30j 

Results: 

In Fig. HKb) the single particle densities of states of 
pure and impurity doped graphene are shown. They 
feature two characteristic singularities around ±3.8eV 
and a V-shaped dip at the Fermi energy. Please note 
that there is a tail of states in the doped system beyond 
the regular band width and also some additional states 
around lo — 0. In Figfljc) the spectral density of all 
plasmon and the spectral density of localized plasmons 
are shown for p = and Uq = 3.7eV. The global den- 
sity of plasmons can be understood in terms of the single 
particle density of states shown in FigQTb). There are 
very few single particle states available around lo = 0. 
Then there are two peaks around lo = ±3.8eU. And, 
at the extrema of the band the single particle density 
of states becomes small again. This implies a quadrat- 
ically increasing spectral density at small lo, a maxi- 
mum around lo = 3.8 — (—3.8) = 7.6eV, and a small 
impurity-dominated contribution around lo = 16eU (see 
Fig. QTc)). Localized modes occur throughout the avail- 
able frequency spectrum. Non-local plasmonic modes are 
most abundant around 8eV, because of the larger phase 
space provided by the van-Hove singularities. The high- 
est energy modes are all localized, i.e. the spectral den- 
sity of local plasmons equals the total spectral density 
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with above 16eV. This is a consequence of the fact that the 
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FIG. 1: (Color online) (a) Single impurity on the graphene 
lattice. Blue corresponds to OeV and red to 3.7eV. The black 
dot indicates the center of the impurity xo. (b) Single particle 
density of states n(u)) = —^Im trQ°(uj) = S(uj — Ep) for 
pure((7o = 0) and doped((7o = 3.7eV) graphene, both with 
H — 0. (c) Spectral density of plasmons in graphene with 
/i = and Uo = 3.7eV. "Global" corresponds to the density 
of states trA(cj) — ~^2 all A a (uS), and "Local" is the spectral 
density of plasmons localized around the impurity trA(ui) = 
^2, , A a (cj), here a runs through localized modes only. 
The arrows indicate the modes shown in Fig[2] 



high energy plasmons are excitations from the lowest en- 
ergy states, which tend to be uniformly spread through- 
out the lattice, to the tail of highest (localized) single- 
particle states, thus generating a very localized charge 
density profile. The high energy values for these localized 
plasmonic modes are set by the very localized nature of 
the impurity problem. Had the impurity potential range 
been much larger, on the scale of lOnra, one would ex- 
pect those localized plasmonic modes with much smaller 
energy, on the scale of few eV. 

In Figs. (Ua)-(d) the spatial profiles of some selected 
localized modes are shown. These are representatives for 
the diversity of localized modes. It is seen that some 
modes have strong dipole characteristics, whereas others 
have a strong quadrupole component. This point high- 
lights the importance of resolving the spatial distribution 
of the induced charge for all possible modes. Previous 
methods would only be capable of detecting selected few 
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FIG. 2: (Color online) (a)-(d) Charge density profiles of some 
localized plasmons. Red (blue) corresponds to negative (pos- 
itive) charges. The plots show the induced charge (An<f) c ')bb 
for the mode a on site b. The frequencies ui and strengths 
A a = (4> a ) T An4> a of each mode are also indicated. 



modes with strong dipole moments.^' 
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FIG. 3: (Color online) Dependence of the spectral density 
trA(u)) of localized plasmon modes on the chemical potential 
fi and the frequency to. The impurity potential is kept fixed 
at Uo = 3.7eV\ The spectral density is shown via the color 
scale. 

The dependence of the intensity of the localized plas- 
mons (for a fixed impurity potential) on the chemical 
potential is shown in Fig. [3] For larger chemical poten- 
tials, more modes are present at lower energies. This is 
a consequence of the occupation of single particle states 
closer to those localized modes at higher energies. When 
the chemical potential has opposite sign with respect to 
the impurity potential, there is some spectral density at 
low energies. These modes stem from the extra single 
particle states around u> = in doped graphene, seen in 
Fig. Hlb). This shows that the chemical potential, which 
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can be experimentally controlled through a gate voltage, 
is an important tuning parameter for achieving certain 
target modes. A limiting aspect of this property though 
is that the frequency of the modes does not change very 
significantly. This happens because the chemical poten- 
tial does not change the single particle states but only 
their occupation. In this sense, the chemical potential 
does not allow one to control the frequency of localized 
plasmons but only their existence. As it will be shown 
below the impurity potential Uq changes the frequency of 
the localized modes more effectively. We emphasize that 
Fig. [3] shows localized plasmonic excitations. 31 1 
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FIG. 4: (Color online) Dependence of the spectral density 
trA(u) for localized plasmons on the impurity potential mag- 
nitude Uo and on the frequency u. The chemical potential is 
kept fixed fj, — 0. The spectral density is shown via the color 
scale. 

In order to change the frequency of modes, one can 
vary the impurity potential, because it actually affects 
the single particle states. Fig. 2] shows that for an ap- 
propriate change of Uq the frequency of some localized 
plasmons can be tuned by about leV. Note also that 
the impurity potential should be at least 2eV in order 
to obtain a significant intensity of localized modes, and 
some modes are strong only within a specific range of the 
impurity potential (see e.g. the features at to ~ 6.0eV 
or w ~ 7.5eV). The impurity potential considerably 
changes the spectral intensity of localized plasmons and 
can therefore be used as a tuning parameter to achieve 
targeted spectral properties, especially when combined 
with variations of the chemical potential. This confirms 
should be expected since it is already known that the im- 
purity changes the single particle density of states. fl9i — 
21| Another striking feature displayed in this figure is the 
symmetry of the intensity of the modes with respect to 
the sign of the impurity potential. This is a consequence 
of the single particle spectrum being symmetric relative 
to fi = 0. Since the plasmon involves electrons and holes, 
its properties do not depend on the sign of the impurity 



potential. 

Conclusions: We have introduced an RPA approach 
which resolves the real space structure of plasmonic 
modes in graphene. This method was used to show that 
impurities induce the formation of nanoscale localized 
plasmonic excitations in graphene sheets. The spatial 
profile, i.e. dipole vs. multipole, of the modes was found 
to vary strongly with the particular resonance. Further- 
more, their frequency and amplitude can be tuned by 
varying the strength of the impurity potential. We also 
studied the effect of varying the chemical potential on 
these modes. It was found that the chemical potential 
can be used to turn them on and off, but it does not affect 
their frequency. This theoretical study is a first step in 
exploring surface enhancement phenomena in graphene 
which may proof useful for nanoscale technologies such 
as molecular sensing. 
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